SOFTWARE GPS BASED INTEGRATED NAVIGATION 



Technical Field of the Invention 

The present invention relates to signal 
5 acquisition and tracking in a GPS system. 



Background of the Invention 

The Global Positioning System (GPS) is a 
satellite-based navigation system developed and operated 

10 by the U.S. Department of Defense. GPS permits land, 
sea, and airborne users to determine their three- 
dimensional positions, velocities, and time. 

GPS uses NAVSTAR (NAVigation Satellite Timing 
and Ranging) satellites. The current constellation 

15 consists of 21 operational satellites and 3 active 

spares. This constellation provides a GPS receiver with 
four to twelve useable satellites from which it can 
receive GPS satellite signals at any time. A minimum of 
four satellites allows the GPS receiver to compute its 

20 GPS position (latitude, longitude, and altitude) and GPS 
system time. Altitude is typically referenced to mean 
sea level. 

The GPS satellite signal from the NAVSTAR GPS 
satellites contains information used to identify the 
25 satellite, as well as to provide position, timing, 
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ranging data, satellite status, and the updated ephemeris 
(orbital parameters) . 

A worldwide system of tracking and monitoring 
stations measures signals from the GPS satellites and 
5 relays information from the signals to a Master Control 
Station. The Master Control Station uses this 
information to compute precise orbital models for the 
entire GPS constellation. This information is then 
formatted into updated navigation messages for each 
10 satellite. 

The GPS receivers contain antennas, reception 
equipment, and processors that are utilized for 
determining position and timing based on satellite 
ranging . 

15 Current designs utilize frequency-specific 

hardware at the front end of GPS receivers to perform 
signal acquisition and tracking functions. Such a GPS 
receiver includes a GPS front end and a system processor. 
The GPS front end receives the incoming analog GPS signal 

20 from a GPS antenna. This analog GPS signal is amplified 
by a low noise amplifier and is filtered by a bandpass 
filter. The low noise amplifier and the bandpass filter 
may be provided in multiple stages, if desired. 
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The analog GPS signal at the output of the 
bandpass filter is downconverted in a frequency 
downconverter by mixing the analog GPS signal with a 
local oscillator signal to shift the carrier frequency of 
5 the incoming analog GPS signal to a lower and more 

manageable frequency band. This downconversion can be 
performed multiple times to bring the frequency of the 
analog GPS signal down in steps to the final desired 
frequency. Each mixing operation produces a high 

10 frequency information band along with the desired lower 
frequency band. Accordingly, each mixing stage of the 
frequency downconverter requires a bandpass filter to 
remove the information in the higher frequency band. A 
frequency reference is provided for the frequency 

15 downconverter. As indicated above, the GPS front end is 
typically implemented in hardware. 

The system processor includes an A/D sampler 
which, in response to the signal from the frequency 
reference, converts the downconverted analog GPS signal 

20 to a digital GPS signal. A GPS processor uses the 
downconverted digital GPS signal to first determine 
ranging and satellite information for the GPS satellites 
in the line of sight of its antenna at the time, and to 
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then determine latitude, longitude, altitude, and/or GPS 
system time. 

The mixing stages and band pass filters 
provided by the frequency downconverter comprise typical 
5 frequency specific hardware in the GPS front end of the 
GPS receiver. The design and development of this mixing 
and band pass filter hardware is an expensive part of the 
GPS receiver. Not only does this mixing and band pass 
filter hardware add recurring cost, but mixing circuit 

10 designs and re-designs add non-recurring cost whenever 
even a small change is made to the frequency 
characteristics of the system. Furthermore, this 
frequency-specific hardware also typically suffers from 
changes as a result of aging. 

15 In order to lower cost and to reduce the 

effects of aging hardware, it may be desirable to remove 
this hardware and replace its function with software. In 
particular, a software solution where signal acquisition 
is executed with frequency domain processing and signal 

20 tracking is executed with time domain processing has 
shown to be a promising technique. Frequency domain 
processing techniques may involve Fast Fourier Transforms 
(FFT) for Wavelet Mult iresolut ion Analysis (WMA) . 
However, implementation of a such software based GPS 
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receiver is challenging. One particular challenge 
encountered in designing a software based GPS receiver is 
in performing the handover between frequency domain 
acquisition and time domain tracking. This challenge 
5 arises because, if software based GPS position 

acquisition is done in the frequency domain with a 
"batch" of data, acquisition can take a non-trivial 
amount of time and, thus, there is a time gap between 
when the data used for acquisition is collected and when 

10 the acquisition solution is actually computed. Moreover, 
implementing signal acquisition in the frequency domain 
may be desirable for reasons other than reducing 
hardware, in which case the time required to hand over 
processing from acquisition to tracking is still a 

15 problem. 

The present invention overcomes one or more of 
these or other problems. 

Summary of the Invention 
20 According to one aspect of the present 

invention, a method performed by a GPS receiver of 
processing GPS data comprises the following: acquiring 
GPS data from at least one GPS satellite, wherein the GPS 
data is useful in determining GPS position of the GPS 
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receiver; performing acquisition functions in the 
frequency domain based on the GPS data, wherein the 
acquisition functions produce acquisition results, and 
wherein the performing of the acquisition functions cause 
5 a delay in handing over the acquisition results to a 

tracking function; propagating the acquisition results 
forward in time by an amount of time substantially equal 
to the delay; and, handing over the propagated 
acquisition results to the tracking function. 

10 According to another aspect of the present 

invention, a method performed by a GPS receiver of 
processing GPS data comprises the following: a) 
acquiring GPS data from at least one GPS satellite; b) 
determining a code offset based on the acquired GPS data; 

15 c) determining a frequency offset based on the acquired 

GPS data; d) determining a change in GPS position of the 
GPS receiver during the determination of the code offset; 
e) determining a change in rate of the GPS receiver 
during the determination of the frequency offset; f) 

20 determining an updated code offset based on the code 

offset determined at b) and the change in GPS position; 
g) determining an updated frequency offset based on the 
frequency offset determined at c) and the change in rate; 
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and, h) handing over the updated code offset and the 
updated frequency offset to a tracking function. 

According to yet another aspect of the present 
invention, a method performed by a GPS receiver of 
5 processing GPS data comprises the following: 

downconverting and sampling a GPS signal from at least 
one GPS satellite to produce digital GPS data; 
determining parameters related to position and rate of 
the GPS receiver from the digital GPS data; determining 

10 a change in the position and rate of the GPS receiver 
between a start of the determination of the parameters 
and an end of the determination of the parameters; 
updating the parameters in accordance with the change in 
the position and rate; and, handing over the updated 

15 parameters to a tracking function. 

Brief Description of the Drawings 

These and other features and advantages will 
become more apparent from a detailed consideration of the 
20 invention when taken in conjunction with the drawings in 
which : 

Figure 1 illustrates a software based GPS 
receiver according to the present invention; and, 
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Figure 2 is a flow chart illustrating exemplary 
software executed by the GPS processor of Figure 1. 

Detailed Description 
5 As shown in Figure 1, a GPS receiver 10 

includes a GPS front end 12 and a system processor 14. 
The GPS front end 12 receives the incoming analog GPS 
signal from a GPS antenna 16. This analog GPS signal is 
amplified by a low noise amplifier 18 and is filtered by 
10 a bandpass filter 20. The low noise amplifier 18 and the 
bandpass filter 20 may be provided in multiple stages, if 
desired . 

The analog GPS signal at the output of the 
bandpass filter 20 is downconverted in a frequency 

15 downconverter 22 by mixing the analog GPS signal with a 

local oscillator signal to shift the carrier frequency of 
the incoming analog GPS signal to a lower and more 
manageable frequency band. This downconversion can be 
performed multiple times to bring the frequency of the 

20 analog GPS signal down in steps to the final desired 
frequency. Each mixing operation produces a high 
frequency information band along with the lower frequency 
band. Therefore, each mixing stage of the frequency 
downconverter 22 may require a bandpass filter to remove 
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the information in the corresponding higher frequency 
band. A frequency reference 24 provides a frequency 
reference for the frequency downconverter 22. 

The system processor 14 includes an A/D sampler 
5 26 which, in response to the signal from the frequency 
reference 24, converts the downconverted analog GPS 
signal to a digital GPS signal. A GPS processor 28 uses 
the downconverted digital GPS signal to first determine 
ranging and satellite information for the GPS satellites 

10 in view of its antenna at the time, and to then determine 
latitude, longitude, altitude, and/or GPS system time. 

Acquisition and tracking can be performed by 
the GPS processor 28 in software. Acquiring and tracking 
the GPS signals principally means that the code and 

15 frequency offsets between satellites and receiver are 
determined and tracked on a regular basis. A typical 
value of how often the tracking loops execute is 1000 
times per second. GPS position is determined from these 
offsets, usually at a rate of 1 Hz, but it can be done 

20 much faster if the processing capability is available in 
the GPS receiver. 

Thus, the GPS processor 28 may be arranged to 
determine the GPS position of the platform supporting the 
GPS receiver 10 from the acquired GPS signals. That is, 
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the GPS processor 28 determines latitude and longitude of 
the GPS receiver/antenna 10/16 from the acquired GPS 
signal. The GPS processor 28 can also determine altitude 
(or depth) if the GPS receiver /antenna 10/16 operates in 
5 three dimensional space. The GPS processor 28 may 
further determine rate and GPS system time from the 
acquired GPS signal. 

The GPS data can be processed by the GPS 
processor 28 in the frequency domain instead of the time 

10 domain. Frequency domain processing techniques may 
involve Fast Fourier Transforms (FFT) for Wavelet 
Multiresolution Analysis (WMA) . Software based signal 
acquisition is faster than hardware based signal 
acquisition, and software processing can eliminate much 

15 of the front end processing. However, while frequency 
domain techniques during signal acquisition are many 
times faster than conventional time domain techniques, 
there still will be a certain amount of time required to 
complete the computation required for acquisition. Thus, 

20 a delay will occur from the end of the batch data 

collection and the completion of the signal acquisition 
phase. This delay causes a problem in handover between 
the acquisition and tracking functions of the GPS 
receiver 10 because one or more of the GPS variables, 
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such as position and/or rate, may very well change during 
this delay. 

The problem of handover between acquisition and 
tracking functions caused by this delay can be 
5 accommodated by using an inertial measurement unit (IMU) 
30 and a navigation processor 32 in the system processor 
14. The inertial measurement unit 30 supplies inertial 
data to the navigation processor 32. Alternatively, an 
inertial navigation system (INS) can be used in place of 

10 the inertial measurement unit 30. As a further 

alternative, both the inertial measurement unit 30 and an 
inertial navigation system can be used to facilitate 
handing over of signal acquisition to tracking. 

Figure 2 is a flow chart illustrating a 

15 software routine 50 that can be executed by the GPS 

processor 28 in order to perform signal acquisition based 
on GPS data received with the GPS signal at the antenna 
16. A block 52 starts the signal acquisition process, 
and a block 54 notes the start time of signal 

20 acquisition. The block 54 also notes the start 

navigation state as supplied by the inertial measurement 
unit 30 and the navigation processor 32. The navigation 
processor 32 computes the navigation state which 
typically consists of the three velocities, three 
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positions, and three angular orientations of the 
platform, such as a vehicle, supporting the GPS receiver 
10, but other values describing the general state of this 
platform can also be included. The three angular 
5 orientations are frequently stored in the form of a 

direction cosine matrix (DCM) typically named C. The DCM 
C is a 3x3 matrix from which the angular orientations can 
be computed. Also, the horizontal positions, typically 
latitude and longitude, are frequently stored in the form 

10 of a DCM named D. The DCM D is a 3x3 matrix from which 
the horizontal positions can be computed. The 
acquisition start time and the current navigation state 
are important for later processing as discussed below. 

Once the acquisition start time and the start 

15 navigation state are noted, a block 56 initiates the 
collection of the batch of data that is required for 
signal acquisition. The block 58 performs signal 
acquisition by determining the code and frequency offsets 
from which the position of the GPS receiver /antenna 10/16 

20 can be determined. The determination of these code and 

frequency offsets in the time domain are well known. The 
determination of these code and frequency offsets in the 
frequency domain are disclosed in "Comparison of 
Acquisition Methods for Software GPS Receiver/' by David 
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M . Lin and James B.Y. Tsui, ION-GPS 2000 conference and 
"Global Positioning System: Theory & Applications , " 
(Volumes One and Two) by Bradford W. Parkinson (Editor) , 
James J. Spilker (Editor), James J. Spilker, Per Enge . 
5 These code and frequency offsets are handed over to the 
GPS signal tracking function as initial values for its 
tracking states. In the present invention, signal 
acquisition is executed in the frequency domain of the 
GPS RF data as indicated by Figure 2. 

10 Because of the time required for signal 

acquisition to be performed by the block 58, there is a 
substantial likelihood that the GPS position and rate of 
the GPS receiver/antenna 10/16 have changed since the 
beginning of the acquisition process performed by the 

15 block 58. Accordingly, a block 60 determines the changes 
in the code and frequency offsets as computed by the 
acquisition block 58. The computation of the changes in 
the code and frequency offsets is done simply by 
computing the changes in range and range rate from the 

20 GPS receiver/antenna 10/16 to GPS satellite between the 
start of acquisition and the completion of acquisition. 
The block 60 uses the start time noted by the block 54 
and the current time to determine the time that has 
elapsed since the start of signal acquisition. The block 
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60 uses the elapsed time as well as the start navigation 
state noted by the block 54 and the end navigation state 
at the end of signal acquisition in order to compute the 
changes in position and rate from the time at which 
5 signal acquisition is completed by the block 58. 

As discussed above, in propagating the position 
and rate solution of the GPS receiver 10 forward in time 
by an amount of time equal to the delay, the block 60 
uses data from the inertial measurement unit 30. The 

10 navigation processor 32 uses the data from the inertial 
measurement unit 30 to calculate the start navigation 
state (position and rate) at the beginning of signal 
acquisition and the end navigation state at the end of 
signal acquisition, and the block 60 uses these start and 

15 end navigation states to propagate the position and rate 
solution for the GPS receiver 10 from the block 58 
forward in time. 

The navigation processor 32 can use classic 
navigation equations to compute the navigation states 

20 (position and rate) . Vehicular angular velocity and 
vehicle non-gravitational acceleration are the inputs 
from the inertial measurement unit 30. 

These equations are given as follows: 
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C = C^o B \ - {p + Q}c 

v = CA B - {2Q + p)v + g 

D = D{p} 



5 where 

C = transformation matrix from the body (or 
platform) to the local vertical frame 

D = transformation matrix from the local- 
vertical frame to the Earth reference frame 
10 v = vehicle velocity relative to the Earth 

g = plumb-bob gravity vector 

Q = Earth' s angular velocity relative to an 

inertial frame 
p = angular velocity of the local-vertical 
15 frame relative to the Earth 

c5 B = vehicle angular velocity relative to an 
inertial frame (in body axes) 

A B = vehicle nongravitational acceleration (in 
body axes) 

20 A = vehicle nongravitational acceleration (in 

local -vertical axes) 
x = skew-symmetric realization of the vector x. 
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It is noted that all vectors (except To and A ) are 
expressed with components in the local-vertical frame. 
As discussed above, the angular orientations are computed 
from the matrix C, and the horizontal positions, 
5 typically latitude and longitude, are computed from the 
matrix D. These computations are well known. 

The navigation processor 32 uses these 
equations to calculate the start . navigation state 
(position and rate) at the beginning of signal 

10 acquisition, and the navigation processor 32 also uses 
these equations to calculate the end navigation state 
(position and rate) at the end of signal acquisition. 
These two navigation states determine the change in the 
NAV solution between the beginning of signal acquisition 

15 and the end of signal acquisition 

The navigation processor 32 may also use other 
aids in addition to the inertial measurement unit 30 to 
compute its navigation solution. Other typical aids 
include radar or barometric altimeter, Doppler radar, and 

20 even GPS. Of course, during the acquisition phase being 
addressed by the present invention, the GPS receiver 10 
will not have any measurements available to provide to 
the navigation processor. 
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Based on this change in NAV solution, a block 
62 computes the updated code offset from the change in 
position of the GPS receiver 10 between the last pass 
through the software routine 50 and the current pass 
5 through the software routine 50, and a block 64 computes 
the updated frequency offset from the change in rate 
between the last pass through the software routine 50 and 
the current pass through the software routine 50. 

The updated code and frequency offsets to be 
10 handed to the tracking function can be computed based on 
the equations below. First, certain quantities are 
defined as follows: 



SVPb: SV position at the beginning of 
15 acquisition as computed with almanac data typically 

computed in Earth Centered Earth Fixed (ECEF) coordinates 

SVPe: SV position at the end of acquisition as 
computed with almanac data typically computed in Earth 
Centered Earth Fixed (ECEF) coordinates 
20 Xb: user position at the beginning of 

acquisition typically computed in Earth Centered Earth 
Fixed (ECEF) coordinates 
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Xe: user position at the end of acquisition 
typically computed in Earth Centered Earth Fixed (ECEF) 
coordinates 

dr: change in range between the GPS 
5 receiver/antenna 10/16 and the satellite during the 

signal acquisition phase = mag(SVPe-Xe) - mag(SVPb - Xb) 

SVVb: SV velocity at the beginning of 
acquisition as computed with almanac data typically 
computed in Earth Centered Earth Fixed (ECEF) coordinates 
10 SVVe : SV velocity at the end of acquisition as 

computed with almanac data typically computed in Earth 
Centered Earth Fixed (ECEF) coordinates 

Vb: user velocity at the beginning of 
acquisition typically computed in Earth Centered Earth 
15 Fixed (ECEF) coordinates 

Ve : user velocity at the end of acquisition 
typically computed in Earth Centered Earth Fixed (ECEF) 
coordinates 

dv: change in range rate between the GPS 
20 receiver/antenna 10/16 and the satellite during the 

signal acquisition phase = mag(SVVe-Ve) - mag (SVVb - Vb) 

co_original : code offset computed by 
acquisition 
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co_updated: the updated code offset value to be 
handed to the tracking function 

fo_original: frequency offset computed by 
acquisition 

5 fo_updated: the updated frequency offset value 

to be handed to the tracking function 

chipping_rate : the number of chips per meter 
that the tracking loop uses for its code generation 

lambda: the wavelength of the GPS signal being 
10 tracked (the carrier frequency wavelength) 

Then the updated code and frequency offsets can be 
calculated as by the following equations: 

15 co_updated = co_original + dr *chipping_rate . 

fo_updated = fo_original + dv/lambda. 

It is noted that the SV position and velocities, which 
are required to compute dr and dv, are normally computed 
20 by the GPS receiver based on almanac data and are 
therefore available . 

Finally, a block 66 hands over the updated 
frequency and code offsets to the tracking function that 
is conventionally performed by GPS receivers. 
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Thus, the outputs of the GPS signal acquisition 
represented by the software routine 50 are updated code 
and frequency offsets for each satellite. The inertial 
data from the inertial measurement unit 30 can be used to 
5 propagate the position and rate solution forward from 
batch data collection to the current time. These 
propagated solutions can then be used to initialize the 
tracking function . 

Certain modifications of the present invention 

10 will occur to those practicing in the art of the present 
invention. For example, as shown above, the GPS front 
end 12 includes the frequency downconverter 22 and such 
bandpass filters as are desirable. However, operations 
of the frequency downconverter 22 and/or bandpass filters 

15 instead may be executed in software by the system 
processor 14 . 

Accordingly, the description of the present 
invention is to be construed as illustrative only and is 
for the purpose of teaching those skilled in the art the 

20 best mode of carrying out the invention. The details may 
be varied substantially without departing from the spirit 
of the invention, and the exclusive use of all 
modifications which are within the scope of the appended 
claims is reserved. 

25 
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